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1. INTRODUCTION 

Due to its physical structure, a riser basically is modeled as a tensioned beam [1], [2], [3], and [4]. In 
[5], an active boundary control that produces a vibration-free for an Euler-Bernoulli beam system was designed. 
Similar use of distributed control can be found in [6]. In [7], the authors used differential evolution optimization 
to search for the best controller model structure and its parameters for beam control problem. The proposed 
controller is able to suppress the beam?s vibration without knowledge of the system. However, the searching 
process is conducted within a set of predefined control structures, no proof of the effectiveness of the control 
was given. 

With efforts to make voltage-source converter (VSC) more efficient in handling distributed parameter 
systems, sliding-mode control (SMC) was given extra flexibility by adding a neural network and fuzzy control 
in [8]. The author yield a control law in the form of a mass-damper-spring system at the boundary of a moving 
string. However, difficulties in selecting proper fuzzy membership functions and a slow convergence speed due 
to online-tuning might be troublesome when applying the aforementioned controls. After accepting that SMC 
is non-analytical in the sliding surface, in the first control structure, a boundary layer was defined that enabled 
fuzzy control by taking a switching function and its derivative as inputs while SMC was activated outside this 
boundary to achieve fast transient responses. A series of papers with applications of SMC to flexible system 
can be found in [9], [10], and [11]. A second attempt was made to design a fuzzy neural network control 
(FNNC) that also employed switching variables as its inputs. The proposed FNNC conducted an online-tuning 
process to regulate fuzzy reasoning to compromise system uncertainties. Both controls resulted in a variation 
of axially moving string tension as the control action. 

In [12], a beam model representing a tensioned riser is investigated, and a boundary controller con- 
sisting of the top-end rise information is designed to achieve exponential stability. Krstic, et al. develop a sys- 
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tematic approach based on backstepping control for beam-type structure in [13] and [14]. In [15], the authors 
proposed a control assisted by a disturbance estimator to guarantee asymptotic stability of an Euler-Bernoulli 
beam system subjected to unknown disturbances. He, et al. In [16], successfully develop a boundary control 
for a flexible riser with vessel dynamics. In [17], the authors introduce control based on Lyapunov’s approach. 
Through Lyapunov’s direct method, the riser’s transverse motion under time-varying distributed loads stability 
is established. A control problem for a coupled nonlinear riser exhibiting longitudinal-transverse couplings is 
investigated in [3]. Analogous applications to flexible systems are evidenced in [18], [2], and [19]. Since the 
surface vessel is always control by a dynamic positioning system in practice [20], [21], [22], [23], [24], and 
[25] the vessel’s motions normally are not considered. The paper deals with the vibration control problem for 
marine risers under environmental disturbances. In addition, the longitudinal-transverse coupling in the riser 
motion in taken into account. Different from [26], the control is formulated without the assumption of positive 
tension. Existence, uniqueness, and convergence of the solutions of the closed-loop system is verified in the 


paper. 


2. MATHEMATICAL FORMULATION 
The riser kinetic energy is specified by 


m -Ti a J (POY) de, (1) 





where u(z,t) is transverse displacements in the X direction and w(z,t) is longitudinal displacement in the 
Z direction. L denote the riser length, mo = pA is the riser oscillating mass per unit length, A is the riser 
cross-section area, and p represents the mass density of the riser. Assuming that the riser is constrained by 
constant tension Po. The riser potential energy is given as 
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where F is the Young’s modulus and Z is the second moment of the riser’s cross section area. The hydrodynamic 
forces can be given as [26] 


Taa) = fume Juls 2st) = oD F ub: 
TaD = -Qı puz(z, t), liane = —OQopv;z(z, t), (3) 


where fuD, fwp and fuL, fwr correspond to the distributed damping and external forces. The work done by 
the hydrodynamic forces acting on the system is calculated as 


L É 
Wy = / fulz, t)u(<, t)dz + J faa Dw tdZ, (4) 
0 0 
The work done by boundary control is 
Wm = Uu (L, thu(L, t) + Uy (L, t)w(L, t), (5) 


where U,,(L,t) and U,,(L,t) are the boundary control forces. The total work done on the system is W = 
Wy + Wm. The extended Hamilton principle is indicated as 


t2 
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For the sake of clear presentation, (z, t) is omitted whenever it is applicable. The kinetic energy variation can 


be written as 
t2 
r é6Tdz = —mMo i [CG J2 E + P Sw) ded (7) 
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where ĝu = ĝv = Ow = O at t = t1, t2 have been used. For the riser under consideration, ball joints arranged 
at both ends (Figure 1) implying that bending free. In addition, the lower end stationed at the well-head. The 
riser dynamics is yielded 
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Figure 1. Riser coordinates 


3. CONTROL DESIGN 
In order to minimize the riser vibration using measured state and applied forces at the top end, we 
consider the following Lyapunov candidate function 
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Since Vt > 0 and u(0, t) = w(0,t) = 0, it can be shown that 
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where yı and ‘yz are positive constants, it can be deduced that 
L p; fF L L o fË 
41? yp; | uaz- | urdz < p | uu,zdz < Lyp | u?dz+ er | ide, (11) 
0 V1 JO 0 0 V1 JO 
L po [E L L p [E 
-414p | w?dz— wedz < pa | wwdz < AL? 2p | w?dz+ wedz. (12) 
0 2 J0 0 0 2 JO 


Vibration attenuation control of ocean marine risers with axial-transverse couplings (Tung Lam Nguyen) 


238 g ISSN: 1693-6930 


The (9) can be lower and upper bounded by 
L L 
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If we select p1, p2, y1, and y2 such that: 
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where c;, for i = 1...4, are strictly positive constants. Differentiating (9) and taking (8) into account yields 
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and noted that — ETu,++(L, t)+Pouz (L, t) +24 u3 (L, t)+EAw,(L, t)u,(L,t) = Uu (L, t) and EAw,(L, t)+ 
Eâ u?(L,t) = Uu (L, t), the boundary controls are designed as follows, 


Uy = —kyu(L,t) = kəu(L, t), Uy = —k3w(L, t) = kawi(L,t), (19) 
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where coefficients k;, for? = 1...4, are strictly positive constants. Substituting the controls (19) into (16) 
gives 
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Remark: It is noted that the authors of [26] use the assumption the riser is alsway stretched in order to conclude 


that fi u2w,dt is positive. This is not the case in practice since the riser can be bulked or stretched according 
to external disturbance. Considering the following term 
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To remove the requirement of positive tension, we use the following property [27] that 
22, 1 
wuz + 7 >0 (24) 


From (20), the designed parameters are selected such that 
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where c1, for t = 5...9, are strictly positive constants. Applying the upper bound of V in (14), (20) can be 
written as 
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Remark: Different from [16], the control design process is carried out in this chapter without any assumptions 
on boundedness of time and spatial derivatives of the riser system. 
Equation (26) can be written as 


where 
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An upper bound of A, can be written as 
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There exists a strictly positive constant € such that the following inequality holds 
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From the lower bound of V, it is shown that 
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Substituting (32)and (33) into (29) gives 
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If € is picked such that ¢ = c — $ is strictly positive, then: 
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Inequality (38) implies that V(t) exponentially converges to nonnegative constant 2Q. Using Inequality A.2 
[26], it can be conclude that all terms |w(z,t)| and |w(z,t)| are bounded and exponentially converge to a 
non-negative constant defined be the value of external disturbances. 
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4. NUMERICAL SIMULATIONS 


At this stage, we illustrate the advantages of the proposed control through a set of simulations. 
The marine riser system parameters are given as in Table 1 The linear current velocity vector in a form of 
V = |[+s, 02 s, 0]* is employed in numerical simulations. The hydrodynamic forces can be given as [26]. 
Simulations are carried out without the proposed control and with the control by set kı = kg = 500. The riser 
displacements in the X and Z directions for uncontrolled and controlled cases are plotted in Figure 2 and Figure 
3, respectively. It can be observed that when the control 1s activated, displacement magnitudes in all directions 
(X and Z) are reduced. The reduction in displacement magnitudes illustrates the effectiveness of the proposed 
control in driving the riser to the vicinity of its equilibrium position. It also can be observed in Figure 4 that 
the control forces required to drive the risers are reasonable for the riser under consideration. 


Table 1. The marine riser parameters 


Nomenclature Description Value 
L Length 1000m 
Do Diameter 0.61m 
D; Diameter 0.575m 
Dyg Diameter 0.87m 
Pu Density 1025kg /m’ 
pon Density 1205kg/m°? 
E Young’s modulus 2 x 10!°kg/m? 
Po Tension 2.15 x 10°N 






RSs 












t(s) : t(s) 
(a) (b) 


Figure 2. The riser’s motions without control: (a) u(z, t) and (b) w(z, t) 


u(m) 











(a) (b) 


Figure 3. The riser’s motions with control: (a) u(z, t) and (b) w(z, t) 
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(a) (b) 
Figure 4. Control input: (a) U„ (L, t), and (b) Uu (L, t) 


5. CONCLUSIONS 

The paper copes with minimizing vibration of the marine riser. After deriving the set of equations 
specifying the riser dynamics, the boundary controller applied at the riser top end is designed thank to Lya- 
punov’s direct method without the assumption of positive tension applied to the riser. The ability in stabilizing 
the riser at its equilibrium position of the boundary control is validated analytically and illustrated numerically. 
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